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ABSTRACT 

The  two-dimensional  HELP  computer  code  has  been  modi- 
fied to  solve  the  interaction  between  buried  explosive  charges 
and  simple  structures  located  above  ground  level.  Three 
calculations  were  performed  and  the  results  compared  favorably 
with  results  from  a concurrent  test  program. 

Dynamic  material  properties  experiments  were  performed 
to  provide  the  necessary  soil  equation  of  state  parameters 
which  are  required  as  input  to  the  numerical  model. 
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I.  INTRODUCTION 


Techniques  for  predicting  the  effects  of  buried  explo- 
sive charges  on  structures  in  contact  with  or  above  the  ground 
surface  are  of  interest  to  engineers  and  designers  working 
in  the  area  of  mine  warfare  technology.  The  computer  code, 
HliLP,^  has  recently  been  modified^  to  solve  two-dimensional, 
axisymmetric  problems  involving  the  detonation  of  2 buried 
explosive  charge,  the  expansion  of  the  detonated  products, 
the  interaction  of  these  gasses  with  the  soil,  the  acceleration 
of  the  soil  and  its  interaction  with  a flat  plate  located 
above  the  ground  and  finally,  the  stress  propagation  through 
the  plate.  The  code  has  been  employed  in  the  theoretical 
portion  of  a joint  theoretical/experimental  investigation  of 
buried  mine  effects  and  the  results  of  three  calculations 
are  presented  herein. 


The  experimental  portion  of  the  investigation  included 
basic  soil  property  tests  to  provide  points  on  the  Hugoniot 
needed  for  determining  input  parameters  for  the  numerical 
model  as  well  as  a test  program  for  the  purpose  of  validating 
the  theoretical  predictions. 

Section  II  of  this  report  discusses  the  experimental 
techniques  employed  to  determine  the  required  soil  properties 
and  presents  the  experimental  matrix  along  with  the  desired 
results.  Section  III  describes  the  test  set  up  and  results 
from  the  buried  explosive  charge  experiments.  In  Section  IV 
the  results  from  three  calculations  are  presented  and  compared 
with  data  from  the  experimental  program.  Concluding  remarks 
are  presented  in  Section  V. 


The  results  reported  herein  indicate  that  significant 
advancement  in  the  state-of-the-art  of  predictive  techniques 
for  evaluating  the  effects  of  buried  explosive  charges  has 
been  attained. 
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II.  SOIL  PROPERTIES  EXPERIMENTS 


'i 


2.1  INTRODUCTION 

Successful  application  of  the  computational  tech- 
niques described  earlier ^ depends  on  the  availability  of 
an  accurate  model  of  the  dynamic  behavior  of  the  soil  under 
which  the  charge  is  buried.  The  particular  equation  of 
state  chosen  is  that  of  Tillotson^-I  which  has  the  form 


P = a + 


Ip  + Ay  + By5 


where  n=p/p  = y + 1,  and  p , a,  b,  I , A,  and  B 
0 0 o 

are  constants.  The  values  of  the  constants  must  be  determined 
either  by  experiment  or  by  analogy  with  other  materials. 

These  parameters  will  be  discussed  below  after  describing 
the  experimental  techniques  and  results. 

Two  soils,  each  in  two  different  density-saturation 

states,  were  studied.  The  first  soil,  McCormick  Ranch  Sand 

from  near  Albuquerque,  New  Mexico,  was  chosen  because  static 

test  data^  and  a limited  amount  of  shock  data  ^ were 

already  available.  The  second  soil,  from  Southwest  Research 

Institute  near  San  Antonio,  Texas  (referred  to  below  as  SwRI 

soil)  was  chosen  because  it  had  been  the  test  medium  for 

previous  experiments  with  buried  explosive  charges.^]  The 

"dry"  state  (p  =1.36  gm/cm3,  w = 8.4*)  was  chosen  as  'he  mean 
o 

state  of  the  wetter  (and  less  dense)  of  two  relatively  dry  groups 
of  soils  used  by  Wenzel  and  Esparza.  ^ The  water  content  and  den- 
sity of  the  "wet"  McCormick  Ranch  sand  were  chosen  to  be  near  those 
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[51  x 

used  in  previous  shock  studies1*  J (p  = 2.12  gm/cm  , w = 12%). 

o 

The  properties  of  the  "wet"  SwRI  soil  were  chosen  to  yield  an  inter- 
mediate density  (1.70  gm/cm^)  at  the  plastic  limit  (19.4%  water). 

2.2  MEASUREMENT  TECHNIQUES 

2.2.1  Static  Properties 

The  plastic  limit  of  the  SwRI  soil  was  determined  using 
ASTM  D424-59  procedures.  The  plastic  limit  of  the  McCormick 
Ranch  sand  was  taken  to  be  14  percent,  the  water  content  used 
by  Peterson  and  Gates  Mazanti  and  Holland  ^ give  a plastic 

limit  of  15  percent  which  is  within  the  accuracy  of  the  technique. 

The  water  content  of  both  materials  as  received  was  deter- 
mined according  to  ASTM  D2216-71  specifications.  The  SwRI  soil 
had  an  initial  water  content  greater  than  8.4  percent  so  it  was 
necessary  to  dry  a portion  of  it  to  mix  with  the  remainder  for 
the  "dry"  samples.  The  water  content  of  the  McCormick  Ranch 
sand  as  received  was  also  8.4  percent.  Each  soil  was  divided 
into  two  batches  and  sufficient  water  added  (if  necessary)  to 
bring  the  water  content  to  the  desired  amount.  Thereafter  the 
soils  were  stored  in  sealed  plastic  bags  except  when  the  soil 
was  being  extracted  for  tests.  The  water  content  was  checked  at 
intervals  and  found  to  be  stable  when  stored  in  this  manner. 

The  desired  initial  density  of  test  samples  was 
achieved  by  pressing  a known  mass  of  soil  into  a known 
volume  of  sample  holder.  For  the  dry  soils  this  could  be  done 
by  hand,  pressing  an  aluminum  anvil  on  the  soil  until  it 
rested  on  the  surface  of  the  ring  containing  the  soil.  For 
the  wet  soils  it  was  necessary  to  use  a hydraulic  press  to 
achieve  the  desired  densities.  Our  past  experience  has  been  that 
this  technique  leads  to  density  variations  of  about  5 percent 
with  higher  densities  at  the  center  of  the  sample  for  soils 
near  their  plastic  limit.  Consequently,  we  used  2 1/2  percent 
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less  soil  than  that  calcul-  r-a  to  give  the  proper  density. 
Thi*.  procedure  gives  ’.lie  desired  density  in  the  center  of 
the  sample  where  the  gages  were  located.  For  dry  soils  the 
variation  from  the  desired  dei  ■■  ity  is  less  than  one  percent. 
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2.2.2  l)ynam;.  . Properties 

The  <;,•  •. -nic  equation  of  state  measurements  below 
40  kbar  were  . '..e  v.  ’n-f  t light  gas  gun  facility.  A 

projectile,  wh*~e  fror>  i is  a flat  plate  of  standard 
material,  is  accelerated  in  the  gun  and,  upon  leaving  the 
barrel,  impacts  the  sample,  sending  a sho^k  wave  into  the 
sample.  Gauges  ai  placed  in  the  sample  to  determine  the 
shock  wave  velocity  and  amplitude.  Analysis  of  the  event 
will  be  described  in  Section  2.2.3  after  describing  the 
experiments  in  more  detail. 

Most  of  the  samples  in  this  study  were  in  targets  of 
two  configurations.  In  the  first  arrangement  (Figure  la)  a 
single  layer  of  the  sample  material  was  placed  between  two 
layers,  one  of  aluminum  and  one  of  plexiglas  (pmma) . At  both 
interfaces  there  was  a manganin  piezoresistive  gauge.  In 
order  to  prevent  individual  soil  particles  from  piercing 
the  gauges  and  changing  their  resistence  a layer  of  epoxy 
impregnated  fiberglas  0.013  cm  thick  covered  tne  gauges. 

As  we  discuss  below,  this  technique  was  not  completely  suc- 
cessful. The  gauges  were  connected  to  a bridge  circuit 
(Pulsar  251A)  to  monitor  their  change  in  resistance  as  the 
shock  wave  passed.  The  pressure  was  determined  from  the 
known  piezoresistive  coefficient  of  manganin,  and  the 
shock  velocity  was  determined  from  the  transit  time  through 
a known  thickness  of  sample.  This  configuration  was  used 
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for  the  initial  experiments  on  SwUI  soil  to  aid  in  designing 
later  shots  and  to  determine  the  approximate  equations  of 
state . 

Most  of  the  gas  gun  experiments  employed  the  target 
configuration  shown  in  Fig.  lb.  In  these  experiments  there 
are  three  piezoresistive  gauges--one  between  the  buffer  and 
the  soil  and  two  more  in  the  soil.  In  one  experiment  a one- 
inch  thick  steel  disc  which  was  used  as  a weight  during  the 
sample  assembly  process,  was  accidentally  glued  to  the  rear 
surface.  Since  the  steel  could  not  be  removed  from  the 
specimen  without  destroying  it  and  since  it  could  not  affect 
the  initial  loading,  the  shot  was  carried  out  as  planned. 

In  this  shot  (number  204)  we  were  able  to  determine  an  approxi 
mate  second- shock  Hugoniot  point.  Unfortunately,  neither  the 
type  of  steel  nor  its  heat  treatment  was  known  so  there  are 
considerable  errors  attached  to  that  datum. 

Our  original  intention  had  been  to  perform  several 
experiments  using  explosively  launched  plates  to  obtain 
equation  of  state  data  in  the  range  from  50  to  150  kbar. 
Because  of  problems  with  gauge  survival  only  one  such  experi- 
ment was  performed.  The  experimental  arrangement  is  shown 
in  Fig.  2.  The  gauges  and  recording  procedures  were  similar 
to  those  in  the  gas  gun  experiments. 


Most  of  the  data  reported  below  were  evaluated  using 
impedance  matching  techniques.  In  cases  where  the  pressure 
calculated  from  impedance  matching  using  the  observed  shock 
velocity  disagreed  markedly  from  that  measured  by  the  manganin 
gauges,  the  former  was  accepted  as  correct.  In  all  such  cases 
the  measured  stresses  were  too  high  and  it  is  thought  that 
this  was  due  to  the  stretching  of  the  gauge  as  it  is  impacted 


Fig.  2- -Experimental  arrangement  for  high  explosive 
equation  of  state  measurements. 
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by  the  individual  grains  in  the  soil.  Such  deformation  is 
a common  problem  in  application  of  piezoresistive  gauges  to 
granular  materials.  A variety  of  encapsulating  media  in- 
cluding kapton,  mylar  and  epoxy  impregnated  fiberglass 
were  investigated  to  try  to  prevent  this.  None  proved  com- 
pletely successful  but  the  best  results  were  obtained  with 
the  fiberglass  packages. 

The  records  from  one  experiment  were  apparently 

unaffected  by  gauge  stretching  in  as  much  as  the  measured  and 

inferred  stresses  were  nearly  identical.  Furthermore,  the 

wave  did  not  attenuate  appreciably  between  the  first  two 

gauges.  These  records  were  analysed  by  the  Lagrangian  gauge 

r 7 | 

analysis  scheme  of  Fowles  and  Williams.  J The  computer  pro- 
gram for  this  analysis  (GAGES)  was  supplied  by  Stanford 
Research  Institute. 

2.3  RESULTS 

2.3.1  McCormick  Ranch  Sand 

A total  of  eleven  gas  gun  experiments  were  performed 
on  McCormick  Ranch  sand.  Two  were  technique  shots  which  pro- 
duced no  useful  data.  Two  failed  due  to  bad  projectile  de- 
sign and  on  a fifth  shot,  scopes  misfired.  One  shot  provided 
both  direct  shock  data  and  reflected  shot  data  so  that,  in 
effect,  there  were  seven  successful  experiments. 

The  available  Hugoniot  data  on  McCormick  Ranch  sand 
are  tabulated  in  Table  I.  The  three  experiments  at  low 
stresses  are  taken  from  Peterson  and  Gates.  ^ 

One  of  the  points  (p  = 32.5  kbars)  in  Table  I is  the 
final  state  from  a Lagrangian  gauge  analysis  of  a pair  of 
stress-time  records  as  mentioned  above. 
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TAB  LI:  I 

MCCORMICK  RANCH  SAND 


0.19 

0.049 

0.32 

1.75 

0.5714 

0 

1.48  1 

i 

0.571 

0.6 

0.072 

0.54 

1.71 

0.5848 

0 

1.48 

- 

0.585 

1.5 

0.060 

1.17 

2.23 

0.4484 

0.12 

2.12 

0.944 

0.379 

7.36 

0.512 

1.06 

2.64 

0.3788 

0.085 

1.36 

a. 843 

0.339 

12.0 

0.617 

1.43 

2.39 

0.4184 

0.085 

1.36 

0.804 

0.386 

26.8 

0.500 

2.53 

2.64 

0.3788 

0.12 

2.12 

0.731 

0.330 

27.4 

0.483 

2.68 

2.59 

0.3861 

0.12 

2.12 

0.29 

0.338 

32.5 

0.620 

2.47 

2.83 

0.3535 

0.12 

2.12 

0.714 

0.303 

33.3 

0.575 

2.75 

2.68 

0.3731 

0.1  2 

2.12 

0.712 

0.326 

(58) 

(0.88) 

(3.1) 

(2.96) 

(0.338) 

0.12 

2.12 

0.660 

(0.293 

P = Hugoniot  pressure  (kbar) 

U •=  Particle  velocity  (mm/ysec) 

D - Shock  velocity  (mm/ysec) 

3 

Pj { = Hugoniot  density  .gm/cm  ) 

Vjj  = Hugoniot  volume  = 1/pj,  (cm3/gm) 
w = Water  content  by  dry  weight  (percent) 

3 

Pq  = Initial  density  (gm/cm  ) 

0 = Specific  volume  cf  water  at  P = Pjj  (cm^/gm) 

2 •? 

Vs  = Specific  volume  of  solids  (cm  /gm) 

Data  in  parentheses  ( ) is  for  a doubly  shocked  state. 
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2.3.2  Southwest  Research  Institute  (SwRI)  Soil 

Six  light  gas  gun  shots  and  one  high  explosive  experi- 
ment were  performed  on  SwRI  soil.  Two  of  the  gas  gun  shots 
yielded  both  Hugoniot  and  reflectel  shock  data,  but  three 
resulted  in  no  usable  data.  These  Jesuits  are  tabulated  in 
Table  II. 

2.4  DISCUSSION 

2.4.1  Experimental  Problems  in  Studying  Shocked  Soils 

Soils  are  extremely  complex  composites  of  \oids,  water, 
and  (usually)  several  solid  phases  (predominately  silicates). 
Because  the  various  romponents  have  markedly  different  im- 
pedcnces,  the  shock  front  in  a soil  will  be  highly  convoluted 
on  the  scale  of  individual  grains.  Furthermore,  just  behind 
the  front  of  the  wave  there  will  be  large  velocity  gradients 
until  the  components  have  reached  mechanical  equilibrium. 

This  feature  makes  it  very  difficult  to  use  in-material  gauges 
in  soils.  Thin  foil  gauges  are  chopped  and  sliced  by  the 
passage  of  a shock  front.  Encapsulating  the  gauges  in  sheet 
material,  such  as  mylar,  kapton  or  fiberglass  improves  their 
performance  but  not  enough  to  permit  the  use  of  pi ezoresistive 
gauges  in  dry  soils.  Even  if  these  materials  prevent  the 
actual  cutting  of  the  gauge  element  by  the  soil  particles  they 
often  will  not  prevent  stretching,  Either  cutting  or  stretch- 
ing of  the  piezoresistive  element  will  give  anomalously  high 
stress  readings  such  as  were  seen  in  many  of  the  shots  re- 
ported here. 

The  use  of  either  ytterbium  piezoresistive  gauges  cr 
magnetic  particle  velocity  gauges  has  been  suggested  to  alle- 
viate this  problem.  Ytterbium,  with  a higher  piezoresistive 
coefficient  than  manganin,  will  give  a larger  signal  at  a given 
stress,  so  that  the  anomalous  signal  due  to  stretching  will  be 
relatively  smaller.  This  is  not  a feasible  solution  at  high 
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TABLE  II 
SwRI  SOIL 


p 

U 

D 

VH 

Vo 

P 

0 

Vo 

Vs 

4.24 

0 24 

1.3 

1.67 

0.5995 

0.085 

1.36 

0.883 

0.575 

7.3 

0.465 

1.15 

2.28 

0.4390 

0.085 

1.36 

0.843 

0.405 

12.6 

0.424 

1.75 

2.24 

0.4456 

0.194 

1.70 

0.800 

0.377 

14. S 

0.467 

1.83 

2.28 

0.4378 

0.194 

1.70 

0.789 

0.370 

(52 

±5) 

(0.87 

±.03) 

(3.52 

±.42) 

(2.26) 

(0.443) 

0.194 

1.70 

(0.669) 

(0.399; 

Symbols  as  in  Table  1.  Data  in  parentheses ()  are  from  high 
explosive  shot  and  may  have  considerable  errors 
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ra t he r than  deform  plastically.  Another,  but  less  promising, 
approach  is  to  try  still  other  encapsulating  media.  Cer- 
tainly the  gauge  can  be  protected  by  very  thick  sheets  on  the 
order  oL  1/16  inch,  but  in  that  case  the  effect  of  the  gauge 
package  on  the  shock  propagation  will  be  great  unless  there 
is  a good  impedance  match  between  the  coil  and  the  gage. 

Perhaps  the  most  promising  solution  is  the  use  of  in- 
material magnetic  velocity  gauges.  Since  these  gauges  can 
have  low  initial  resistance  without  losing  signal  level, 
they  can  be  made  of  thick  enough  material  to  withstand  the 
velocity  gradients  in  soils.  The  main  difficulty  with  these 
gauges  is  the  requirement  that  moving  conductors  must  remain 
remote  from  the  target  and  this  limits  the  impedence  of  fryers 
that  can  be  used. 

2.4.2  Calculation  of  Poreless  McCormick  Rar.ch  ^and  Properties 

In  order  to  interpret  the  above  data,  obtained  from  dif- 
ferent soils  and  water  contents,  in  a unified  manner  let  us  use 
a simple  mixture  theory  to  extract  the  "Hugoniot"  of  the  solids 
from  the  soil  results.  To  do  this  we  assume  complete  pore 
collapse  and  assume  that  the  water  (if  present)  and  the  solids 
both  attain  the  specific  volume  appropriate  to  the  separate  phases 
at  the  pressure  attained  in  the  soil. 

Symbolically, 

V„(P)  ■ (vsCP)  * Mw  VK(P))/(l  * mJ  (2) 

where  Vjj (P) , VS(P),  and  VW(P)  arc  the  specific  volume  of 
the  so. i, solids,  and  water,  respectively,  shocked  to  pres- 
sure P and  Mw  is  the  water  content  in  percent  of  dry 
weight.  The  Hugoniot  of  water,  shown  in  Figure  3b,  is  a combi- 
nation of  data  by  Lysnc^  and  the  tabulation  of  Riney,  et  al.  ^ 
Then 

V’>  ■ (l  ♦ Mw)  Vh(P)  ■ Mw  VK(P)  (3) 
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We  have  done  this  calculation  for  all  the  data  avail- 
able on  these  two  soils  and  the  resulting  volume  of  solids 
is  shown  as  a function  of  Hugoniot  pressure  in  Fig.  3. 

All  of  the  McCormick  Ranch  sand  data  above  1 kbar  lie 
near  the  line  ( P)  = 0.382  - 0.0019P  where  V has  the  units  of 
cc/gm  and  P is  in  kbars.  We  take  this  to  indicate  that 
closing  of  pores  in  this  material  takes  place  at  very  low 
stresses.  By  using  this  expression  for  V (P)  we  show  below 
that  we  can  fit  most  of  the  data  well  in  the  shock  velocity- 
particle  velocity  plane. 

The  data  for  the  SwRI  soil  all  lie  to  the  right  of  the 
line  used  to  represent  McCormick  Ranch  sand  as  seen  in  Figure 
3.  For  t lie  wet  soil  the  deviation  is  within  the  scatter  of 
the  data  for  the  McCormick  Ranch  sand  alone.  However,  if  we 
consider  the  difference  to  be  real  then  two  interpretations 
are  possible.  The  Hugoniot  of  the  solids  may  be  the  same  as 
in  McCormick  Ranch  sand,  but  the  SwRI  soil  may  resist  complete 
pore  collapse  to  fairly  high  pressures.  Alternately  the  solids 
in  the  SwRI  soil  may  be  stiffer.  We  show  the  Hugoniot  for  a 
pure  quartz  rock^^as  an  example.  Even  in  this  case  we  must 
conclude  that  the  dry  SwRI  soil  does  not  collapse  until  at 
least  7-8  kbar. 


2.4.3  Implication  of  Poreless  McCormick  Ranch  Sand 

I f we  use  liq.  (2)  and  the  curves  in  Figs.  3a  and  3b  we 
can  calculate  the  expected  Hugoniots  for  the  void- free  water 
bearing  soils  shown  in  Fig.  4.  Solid  curves  are  for  the 
McCormick  Ranch  sand  solids.  Curves  are  shown  for  all  of  the 
water  contents  shown  in  Tables  I and  II  and  each  data  point  is 
connected  to  the  appropriate  curve  by  a straight  line.  A 
comparison  of  the  calculated  void- free  Hugoniots  with  experi- 
mental data  allows  one  to  estimate  the  degree  of  pore  collapse 
that  has  occurred  in  the  soil. 


13 


uisamdte&aia&tiatiM 


■"  ^vmMwrST’svcwi'v, 


The  same  data  as  in  Fig.  4 are  shown  in.  particle  velocity- 
shock  velocity  space  in  Fig.  5.  In  this  figure  we  have  also 
indicated  the  zero  pressure  sound  speed  for  McCormich  Ranch  sand 
calculated  from  the  hydrostatic  compression  curve. ^ There 
is  excel!  mt  agreement  between  the  observations  and  the  curve 
shown.  Also  shown  is  the  linear  D-U  relation  which  has  been 
used  for  the  NTS  Playa  soil  discussed  in  Subsection  2.4.4.  As 
can  be  seen,  it  is  in  good  agreement  with  our  30  kbar  data  on 

McCormick ‘Ranch  sand  but  agreement  is  poor  at  lower  pressures. 

The  very  high  compressibility  shown  by  the  solids  in 

McCormick  Ranch  sand  up  to  35  kbars  cannot  be  used  to  estimate 

V.(P)  for  very  high  pressures.  The  scatter  of  the  data  in 
Fig.  4 is  sufficient  to  preclude  an  estimate  of  the  curvature. 
Nonetheless,  we  can  be  certain  that  before  stresses  of  100  kbar 
are  reached  the  curvature  would  become  evident.  Because  of 
the  high  initial  compressibility  we  should  expect  the  effect 
to  be  fairly  pronounced. 

The  very  high  initial  compressibility  of  the  solids  in 

McCormick  Ranch  sand  may  indicate  that  there  is  more  water  in 

the  samples  than  supposed.  Inter- lamellar  water  in  clay 
minerals  might  account  for  such  an  effect.  This  water  is 
often  difficult  to  remove  completely  at  moderate  tempera- 
tures. We  would  expect  clean  sands  to  be  stiffer  as,  indeed, 
they  are. 

Using  the  curves  of  Fig.  3 we  can  derive  the  parameters 

A and  B of  the  Tillotson  equation.  The  first  term  is  not 

important  at  the  stress  levels  studied  here.  Values  of  a,  b 

2 

and  I have  been  chosen  so  that  the  terms  a + b / (I/IQn  + 1) 

equals  a typical  Gruneisen  parameter  for  soils.  The  values  are 

given  in  Table  III. 


2.4.4  Revised  Equation  of  State  of  NTS  Playa  Soil 

At  very  high  stresses  (in  excess  of  100  kbar)  we  expect 
that  the  silicates  in  the  soil  will  transform  to  new  struc- 
tures with  silicon  in  six  coordination  with  oxygen  instead 
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of  the  normal  four  coordination.  Such  structures  have  been 
observed  in  most  silicates  that  have  been  studied. In 
light  of  this  possibility  we  have  reassessed  the  equation  of 
state  of  NTS  playa  soil  used  in  our  earlier  reports.  The 
revised  equation  of  state  is  shown  schematically  in  Fig.  6. 

We  have  assumed  that  a phase  change  occurs  between  ISO  kbar 
and  250  kbar.  The  result  is  a much  closer  fit  to  the  experi- 
mental data  that  also  agrees  with  predictions  based  on  crystal 
chemical  considerations  and  previous  experience  with  silicates, 
Although  we  have  not  done  so  here,  it  would  be  possible  to 
use  hq.  (3;  to  determine  the  equation  of  state  of  the  solids 
in  the  NTS  playa  soil  to  determine  how  they  compare  with 
available  data  on  silicates. 


2.5 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  investigation  of  the  dynamic  properties  of  two 
soils  lias  led  to  the  following  conclusions: 

• A simple  mixture  theory  can  be  applied  consistently 
to  the  equation  of  state  of  McCormick  Ranch  sand 

in  three  states  of  porosity  and  saturation,  two 
from  the  j -esent  study  and  one  from  ref.  5. 

• The  sol  ids  in  McCormick  Ranch  sand  show  an  apparent 
initial  density  of  2.62  gm/cm3  and  an  apparent 
initial  bulk  modulus  or  20u  kbar. 

• In  incompletely  saturated  McCormick  Raich  sand 
(p  = 1.36  gm/cm3,  w = 8.5  percent  dry  weight) 
collapse  of  pore  volume  occurs  at  a stress  less 
than  7 kbar. 

• Completely  dry  McCormick  Ranch  sand  has  open  pore 
space  at  pressures  at  least  as  high  as  0.6  kbar. 

• Within  experimental  uncertainty  Southwest  Research 
Institute  (SwRI)  soil  has  an  equation  of  state 
consistent  with  the  same  solids  as  inferred  for 
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McCormick  Ranch  sand,  although  the  limited  data 
match  better  the  behavior  of  pure  quartz. 

• In  incompletely  saturated  SwRI  soil  (p  * 1.36  gm/cm3, 

o 

w = 8.5  percent  dry  weight)  complete  collapse  of  pore 
volume  occurs  at  stresses  between  4 and  7 kbar. 

In  order  to  achieve  a more  complete  understanding  of 
the  behavior  of  the  properties  of  soils  in  the  stress  range 
of  interest  to  explosive-soil  interactions  we  make  the 
following  recommendations: 

• A better  gauge  technique  than  used  here  is  required 
for  equation  of  state  studies  especially  in  dry 
soils.  In-material  magnetic  particle  velocity  gauges 
offer  the  best  hope  for  success. 

• Measurements  on  dense  (p  =2.12  gm/cm3)  McCormick 

o 

Ranch  sand  at  about  15  kbar  and  between  50  and  100 
kbar  are  needed  to  determine  the  curvature  of  the 
pressure  volume  relationship  of  the  solids. 

• Measurements  of  the  equation  of  state  of  SwRI  soil 
at  high  stresses  would  be  helpful  in  determining 
the  equivalence  (or  lack  thereof)  between  SwRI  soil 
and  McCormick  Ranch  sand. 

• Measurements  of  the  equations  of  state  of  both  soils 
at  stresses  below  7 kbar  would  be  requiied  to  deter- 
mine the  details  of  pore  collapse.  To  date  no 
theoretical  model  is  available  to  predict  this  be- 
havior in  soils.  Such  a model  would  aid  immensely 
the  understanding  of  low  pressure  soil  behavior. 
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III.  T1:ST  PROGRAM 

This  work  is  an  outgrowth  of  an  earlier  investigation 
of  the  effects  of  buried  high  explosive  charges  on  steel  plates 
conducted  for  MERDC  by  the  Southwest  Research  Institute 
(SwRI)  (Wenzel  and  Esparza^).  In  order  to  improve  the 
capability  to  predict  effects  of  land  mines  on  vehicles, 
a computer  code  simulating  the  explosive-soil-target  interactions 
has  been  developed.  As  part  of  this  latter  effort,  a series 
of  twelve  model  experiments  were  performed  using  precisely 
characterized  soils  to  provide  data  as  a check  on  the  computed 
results.  There  were  two  main  reasons  for  this  approach: 

(1)  there  was  considerable  variation  in  the  properties  of 
the  soil  used  by  SwRI,  and  (2)  there  was  no  dynamic  equation- 
of- state  data  for  the  SwRI  soil  so  that  the  computer  model 
used  McCormick  Ranch  Sand,  an  entirely  different  soil. 

McCormick  Ranch  Sand  was  used  in  all  of  the  present  model 
experiments. 

3.1  BACKGROUND 

The  SwRI  report  has  been  a major  influence  in  this 
work  and,  consequently,  the  primary  thrust  of  this  background 
section  will  be  directed  toward  that  report.  The  brief 
summary  of  earlier  work  is  based  on  the  discussion  of  Wenzel 
and  lisparza!^ 

3.1.1  Experiments  with  Unburied  Explosive  Charges 

The  earliest  systematic  investigations  of  the  effects 
of  explosive  charges  on  vehicles  are  reported  in  several 
BRL  reports  (Hoffman  and  Mills^^;  Goodman^^;  Jack  and 
Armendl ) . These  investigators  measured  pressure  and 
impulse  imparted  to  metal  plates  separated  from  explosive 
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charges  by  air  or  partial  vacuum.  These  data  were  used  by 


Cockrell,  et  al.  and  Wenzel,  et  al. 


to  estimate 


survivability  and  vulnerability  of  tanks  to  land  mines. 


The  loading  of  a metal  plate  from  an  explosive  charge 
in  air  is  accomplished  by  two  means.  First,  a large  ampli- 
tude, short-duration  pulse  is  produced  by  a shock  wave 
travelling  through  the  air.  This  initial  loading  contributes 
very  little  impulse  to  the  plate  because  of  its  short  duration. 
The  main  portion  of  the  total  impulse  to  the  plate  is  delivered 
by  the  expanding  detonation  products  in  a somewhat  lower 
pressure  but  much  longer  duration  pressure  pulse.  Jack  and 
Armendt^6^  have  observed  that  the  initial  sharp  peak  is 


Armondt 1 J have  observed  that  the  initial  sharp  peak  is 
virtually  absent  in  experiments  at  partial  vacuum. 


3.1.2  Experiments  with  Buried  Explosive  Charges 


The  effects  of  burial  on  blast  effects  was  first 


investigated  by  Kincheloe.1 


Comparing  his  results  with 


Goodman' s , he  concluded  that  buried  charges  are  more 
efficient  than  charges  in  air,  and  that  the  impulse  delivered 
is  dependent  on  charge  mass,  separation  distance  and  depth 
of  burial. 


Wenzel  and  Esparza1  i conducted  an  extensive  and 
systematic  investigation  of  the  "pressures  and  impulses  at 
close  distances  from  explosive  charges,  buried  and  ii;  air." 
They  investigated  both  spherical  charges  and  pancake  charges 
with  a diameter  to  thickness  ratio  of  3.16.  They  also 
conducted  a small  number  of  experiments  with  cylindrical 
charges  in  air.  Their  charges  varied  from  0.057  kg  (0.125  lb) 
to  0.68  kg  (1.5  lb),  and  they  used  standoff  distances  from 
0.12  ni/kg1/3  (0.3  ft/lb1/3)  to  0.40  m/kg1/3  (1.0  ft/lbl//3). 


» 


The  SwRI  investigators  measured  impulse  by  the  BRL  plug 
technique.  Pressures  were  measured  by  a Hopkins  on  bar  technique 
in  which  high  strength  steel  rods  [nominal  yield  at  24  kbar 
(350  ksi)]  were  inserted  into  holes  drilled  through  thick  steel 
plates.  The  rods  were  separated  from  the  plates  by  sleeves 
of  silicone  rubber.  The  pressure  history  of  the  end  of  a rod 
was  transmitted  up  the  rod  and  measured  by  a strain  gauge  on 
the  surface  of  the  rod.  The  stress-strain  relation  of  the  rod 
then  permits  conversion  of  strain  to  stress.  This  relation  was 
determined  both  statically  and  dynamically,  although  it  is  not 
clear  from  their  report  if  the  calibration  extended  over  the 
entire  range  of  observations. 


Briefly  summarized,  their  results  on  buried  explosives  were: 


• Peak  stress  falls  off  more  rapidly  away  from  the 
axis  of  the  experiment  for  pancake  charges  than 
for  spherical  ones. 

• There  was  little,  if  any,  observed  effect  of 
scaled  depth  of  burial  for  scaled  depths  of 

2 and  4 inches  (5  and  10  cm).  For  a scaled 
depth  of  10  inches  (25  cm)  there  was  an  increase 
in  peak  stress. 

• Axial  pressures  produced  by  pancake  charges 
averaged  2.5  times  those  produced  by  spheres. 

• Buried  spherical  charges  delivered  considerably 
greater  peak  stress  than  corresponding  charge. - 
in  air.  However,  for  the  pancake  charges  there 
was  little  apparent  difference. 

• Water  content  of  the  soil  has  a significant 
effect  on  the  pressure  delivered  to  the  plate, 
with  wetter  soils  producing  higher  pressures. 


ii 
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3.1.3  Systematics  in  the  SwRI  Data 

In  the  initial  stages  of  the  present  work  we  conducted 
a detailed  re-examination  of  the  SwRI  report  just  described. 

The  goal  of  this  effort  was  to  identify  any  weaknesses  in  their 
techniques  in  order  to  improve  the  quality  of  our  own  experi- 
ments. In  the  course  of  this  re-examination  we  have  discovered 
some  systematic  effects  apparently  overlooked  in  that  report. 

Wenzel  and  Lsparza  tabulate  the  soil  water  content  and  bulk 

density  for  all  their  buried  experiments  except  for  two  shots 

in  "mud."  These  data  are  plotted  in  Figure  7 . In  a plot  of 

this  nature  it  is  readily  apparent  that  the  soils  used  fall 

predominantly  into  one  of  two  groups.  The  larger  of  these 

two  groups  has  water  content  (w)  between  5.14  and  11%  of  dry 

3 3 

soil  weignt  and  density  (p)  from  1.18  gm/cm  to  1.51  gm/cm 

3 

with  a mean  of  w=  8.4%  and  p=  1.36  gm/cm  . A smaller,  more 
compact  group  centered  at  w * 3.4%  and  p=  1.53  gm/cm3  has 
2.5%<  w<  4.5%  and  2.48  gm/cm3<p<  2.62  gm/cm3.  The  large  scatter 
of  soil  states  is  undoubtedly  due  to  lack  of  control  of  water 
content  and  density  of  the  soil  directly.  The  SwRI  workers 
compacted  all  soils  by  a similar  process  but  soils  of  different 
water  contents  will  respond  to  that  process  differently. 


The  wide  variability  of  soil  conditions  studied  by  SwRI 
permits  investigation  of  the  effects  of  soil  properties  on 
blast  effects.  Let  us  restrict  our  discussion  to  soils  falling 
into  cither  of  the  two  groups  described  above  and  reject  other 
soils  as  "anomalous."  If  we  then  plot  the  ratio  of  observed 
peak  stress,  , to  the  peak  stress  for  the  same  charge  at 


the  same  distance  in  air,  P , as  a function  of  soil  density 

a 


we  find  a definite  positive  correlation  for  spherical  charges 
(Figure  8a)  but  no  apparent  correlation  for  pancake  charges 
(Figure  8b).  So  we  may  ascribe  some  of  the  scatter  of  the 
SwRI  data  to  variations  of  density. 
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uter  content  plotted  versus  density  for  the  soils 
sed  in  the  experiments  reported  by  Wenzel  and 
sparzal6!.  "Wet"  and  "dry"  groups  are  encircled. 
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Next,  we  note  that  for  nine  of  their  soil  conditions 
Wenzel  and  Esparza  determined  the  sound  speed.  Four  of  these 
cases  were  for  spherical  charges  where  peak  stress  on  axis  was 
also  measured.  In  these  cases  we  can  calculate  the  acoustic 
impedance  ZQ  = P0cQ.  If  the  normalized  peak  pressure  (on  axis) 
is  plotted  against  the  impedance  we  see  an  even  stronger  posi- 
tive correlation  (Figure  9),  It  seems  likely,  therefore,  that 
much  of  the  scatter  in  Figures  8a  and  8b  may  be  due  to  varia- 
tions in  the  sound  speed  of  the  soil.  Based  on  the  scaling 
relations  of  Wenzel  and  Esparza,  this  result  is  not  at  all 
unexpected. 


It  is  interesting  to  note  that  the  actual  depth  of 
burial  seems  to  ^ more  closely  related  to  measured  peak 
stress  than  the  scaled  depth  of  burial.  Figure  10a  is  a plot 
of  real  depth  of  burial  versus  the  normalized  pressure  for 
spherical  charges.  There  is  a weak,  but  nonetheless  distin- 
guishable, effect  of  real  depths  of  burial  on  the  peak  pressure 
as  indicated  by  the  dashed  line  in  the  figure.  When  scaled  depth 
of  burial  is  used  as  the  abscissa  as  in  Figure  10b  no  systematic 
relation  can  be  seen.  We  conclude,  therefore,  that  peak  pressure 
is  more  closoiy  related  to  real  depth  of  burial  than  to  scaled 
depth  of  burial  for  soils. 
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figure  9.  Normalized  peak  pressure  plotted  versus 
impedance  for  various  soil  conditions 
used  in  the  experiments  conducted  by 
Wenzel  and  Esparzat^]. 
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Figure  10a.  Plot  of  real  depth  of  burial  versus  normalized 
pressure  for  spherical  charges.  Dashed  line 
indicates  trend  of  data. 


Figure  10b.  Plot  of  scaled  depth  of  burial  versus  normalized 
pressure  for  spherical  charges. 
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3.2  APPROACH 


3.2.1  General  Experiment  Plan 


In  the  current  effort  twelve  experiments  have  been  con- 
ducted with  buried  high  explosive  charges--two  charge  confi- 
gurations were  used,  spherical  and  3:1  (diameter:  thickness) 
pancake,  and  two  different  soil  conditions  were  used.  Three 
experiments  were  done  for  each  combination  of  charge  and  soil. 

In  all  of  the  experiments  the  charge  was  buried  under  2.54  cm 
(1.00  inch)  of  the  test  soil.  Stress  histories  were  measured 
in  a mild  steel  plate  suspended  over  the  soil  so  its  bottom 
surface  was  at  a scaled  height  of  0.13  m/kg^^  (0.3  ft/lb^^). 

A cross-section  of  the  test  set-up  is  shown  in  Figure  11  for  both 
configurations. 


The  explosive  used  in  these  experiments  was  C-4.  The 
initial  plan  had  been  to  use  1,15  kg  (2-1/2  pounds)  for  each 
charge  whether  spherical  or  pancake.  However,  a 10.3  cm  (4-1/16 
inch)  diameter  spherical  mold  which  held  0.95  kg  (2.1  lb) 
of  C-4  was  already  available,  so  it  was  decided  to  take  advantage 
of  scaling  relations  and  use  the  smaller  charge. 

3.2.2  Soil  Preparation 

The  main  purpose  of  these  experiments  was  to  provide  ex- 
periment:! data  with  which  to  check  the  computer  modelling  of 
buried  charges  described  previously.  Consequently,  it  was  of 
utmost  importance  that  the  soil  used  be  well  characterized. 

The  soil  chosen  was  McCormick  Ranch  sand  from  near  Albuquerque, 
New  Mexico.  This  soil  has  been  investigated  in  both  static 
tests  (Mazanti  and  Holland  ^)  and  dynamic  tests  (Peterson 
and  Gates  ^).  Since  only  a limited  amount  of  soil  was 
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Figure  lib.  Cross-section  of  experimental  arrangement 
for  buried  erplosive  discs. 
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available  to  us  (Stanford  Research  Institute  (SRI)  supplied 
us  with  about  40  kg  of  soil),  it  was  decided  to  use  the 
test  soil  only  above  and  around  the  spherical  charges  and  only 
above  the  pancake  charges.  This  decision  involved  the  trade- 
off between  having  a homogeneous  soil  throughout  the  test  bed 
and  having  a very  precisely  known  soil  where  it  should  have  the 
greatest  effect.  Clearl  , we  have  considered  the  latter  to 
be  the  more  important. 


The  natural  soil  at  the  Green  Farm  Test  Site  is  very  rocky, 
so  enough  soil  was  sifted  through  a 3-mm  (1/8-inch)  screen  to 
conduct  our  tests.  Since  this  rather  tedious  procedure  haci  to 
be  resorted  to,  it  was  decided  to  perform  the  entire  test  in 
a "sand  box"  0.9  in  x 0.9  m x 0.3  m (3  feet  square  by  1 foot 
high).  These  dimensions  are  sufficient  that  edge  effects 
were  unimportant  during  the  period  of  taking  data.  This  pro- 
cedure considerably  shortened  the  process  of  setting  up  the 
experiment . 


The  two  soil  conditions  studied  were  density  p = 1.36  gm/ 


cin  and  water  content  w = 8.5%  of  dry  weight,  and  p = 2.12  gm/ 
3 


cm  and  w = 12%  of  dry  weight.  The  former  is  already  familiar 
to  tiie  reader  as  the  mean  of  the  larger  of  the  two  groups  of 
soil  conditions  studied  by  Wenzel  and  Esparza.^  It  is  also 
rather  close  to  the  in  situ  water  content  of  the  local  soil 


in  late  summer  (7.5%  of  dry  weight).  We  chose  this  soil 
condition  with  the  hope  that  our  results  could  be  compared 
directly  to  those  of  the  earlier  study.  The  latter  soil 


condition  is  that  used  by  Peterson  and  Gates ^ in  their 


low  velocity  gas  gun  experiments.  It  is  near  the  state  studied 
by  Mazanti  and  Holland^ 
w = 11.4%  of  dry  weight. 


r 4 1 3 

by  Mazanti  and  Holland1  1 who  used  p = 1.87  gm/cm  and 
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Soil  was  received  in  three  batches,  two  for  the  final 
experiment  only.  The  first  batch  as  received  from  SRI  had 
a water  content  of  8.41.  The  entire  volume  of  soil  was  sifted 
through  a 1-mm  siev  to  eliminate  the  few  small  pebbles  present. 
This  material  was  us<_  , for  both  the  equation-of-state  studies 
described  previously  a*  d the  present  experiments.  For  the 
drier  soil  condition  this  sifted  soil  was  used  directly,  whereas 
water  was  added  for  the  wetter  condition  to  bring  the  water 
content  to  12% . An  appropriate  mass  of  soil  was  then  loaded 
into  a mold  of  a shape  to  conform  with  the  charge,  and  compressed 
to  the  proper  density.  The  wet  soil,  which  was  cohesive,  was 
then  removed  from  the  mold  and  placed  in  a plastic  bag  to  pre- 
vent water  loss.  The  dry  soil,  which  was  not  cohesive,  was 
retained  in  the  mold  until  buried  when  the  surrounding  soil 
would  support  it. 


3.2.3  Measurement  Techniques 


In  all  of  the  experiments,  stress  histories  were  measured 
at  various  locations  within  the  steel  plate  using  manganin 
piezoresistive  gauges.  Initially,  an  attempt  was  made  to  measure 
the  velocity  of  the  upper  surface  of  the  plate  by  a Moire  tech- 
nique. These  attempts  were  not  successful  because  of  the  low 
velocity  of  the  plate  but  are  summarized  below,  since  the  tech- 
nique would  be  useful  for  situations  for  which  the  rear  free 
surface  velocity  of  the  plate  were  higher.  In  the  final  ten 
experiments  the  surface  velocity  was  measured  by  shorting  pins 
above  the  top  of  the  plate  over  the  center  of  the  charge. 


3. 2. 3.1  Manganin  Gauges 


Stresses  within  the  steel  plates  were  measured  by  manganin 
gauges.  This  technique  was  chosen  for  several  reasons.  First, 
because  of  the  high  density  of  the  soil  in  some  of  our  experiments 
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we  expected  peak  stresses  to  exceed  the  yield  strength  of  any 
available  steel.  This  eliminated  the  use  of  a Hopkinson’s 
bar  technique  such  as  that  used  by  Wenzel  and  Esparza  (1972) . 
Second,  our  personnel  have  considerable  experience  in  both 
the  use  and  calibration  of  these  gauges.  Finally,  even  though 
there  are  some  non-linearities  at  low  stresses  (below  30  kbar) 
aud  considerable  hysteresis  on  unloading,  both  these  conditions 
have  been  sufficiently  studied  so  as  not  to  hinder  interpretation 
of  the  records. 


In  order  to  emplace  the  gauges  within  the  steel  plates 
each  plate  was  made  up  of  a stack  of  two  or  more  plates  with 
gauges  between  the  plates.  The  gauges  were  first  glued  to  the 
lower  plate  with  C-7  epoxy  and,  after  that  epoxy  had  cured, 
the  entire  joint  between  plates  was  filled  with  C-7.  This 
epoxy  layer  was  about  0.4  mm  (0.015  in.)  thick.  Two  gauge 
layouts  were  used.  In  the  majority  of  experiments  a two- layer 
plate  was  used  with  all  the  gauges  located  as  in  Figure  12  in 
the  same  plane.  We  refer  to  this  layout  later  as  H (horizontal). 
In  one  of  these  experiments  (2144)  the  lower  plate  was  5 cm 
thick.  In  three  experiments  (one  for  each  configuration  except 
a spherical  charge  buried  in  dry  soil)  the  plate  consisted  of 
five  plates  with  a gauge  at  each  joint  on  the  axis  of  the 
charge  (sec  Figure  13).  These  experiments  were  conducted  to 
determine  the  amount  of  attenuation  of  the  peak  stress  as 
the  pulse  traveled  through  the  steel  plate.  We  refer  to  this 
' atcr  as  the  V (vertical)  layout. 


The  gauges  used  were  Pulsar  FM50A  gauges.  These  gauges 
are  mounted  between  fiberglass  and  kapton.  The  gauge  itself 
is  only  about  0.002  cm  thick  but  the  entire  package  thickness  is 
controlled  by  internal  solder  joints  and  is  about  0.4  mm  (15 

mils) thick.  The  gauges  were  used  in  a two-lead  arrangement 
and  were  powered  by  a Pulsar  Model  251A  which  also  contains 
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the  bridge  circuit  to  convert  changes  in  gauge  resistance  to 
changes  in  voltage.  Prior  to  each  shot  the  circuits  were  cali- 
brated by  putting  a known  resistance  in  series  with  the  gauge. 

3. 2. 3. 2 Moir£  Fringe  Measurements  of  Surface  Velocity 

Use  of  Moire  fringes  is  a standard  technique  for  the  measure- 
ment of  strain  in  "static"  tests.  We  give  a very  brief  discus- 
sion of  the  general  technique  here  and  refer  the  reader  interested 
in  details  to  the  literature  (e.g.,  Dove  and  Adams^0^).  When 
two  grids  are  superimposed,  a new  pattern  of  alternating  light 
and  dark  bands  can  be  seen.  This  effect,  known  as  the  Moire 
effect,  is  caused  by  overlapping  of  positive  and  negative  features 
of  the  two  grids.  If  one  of  the  grids  is  deformed  the  "fringes" 
show  the  deformation.  It  is  the  mechanical  analog  of  optical 
interference  phenomena.  Such  an  effect  can  also  be  produced 
by  superposition  of  a screen  and  its  shadow. 


We  had  hoped  to  use  this  technique  to  measure  the  motion 
of  the  upper  surface  of  the  plate  in  our  experiments.  Initial 
feasibility  experiments  showed  that  the  method  was  a good  one 
in  principle.  A 1/2- in. -thick , 12-in. -square  steel  plate  was 
placed  on  top  of  a 1/2-pound  sphere  of  C-4.  The  upper  surface 
of  the  plate  had  a grid  of  1/8-in.  lines  on  1/4-in.  centers. 
This  setup  was  viewed  with  a high-speed  framing  camera  through 
a similar  grid.  The  resulting  pictures,  obtained  when  the 
charge  was  detonated,  show  very  clearly  the  deformation  of  the 
steel  plate.  The  results  of  this  experiment  are  given  in 
Figure  14.  Note  that  the  surface  velocity  in  this  case  is  about 
170  m/sec. 


In  our  first  two  buried  charge  experiments,  we  attempted 
to  use  this  same  technique  to  measure  the  surface  velocity. 
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The  optical  arrangement  is  shown  in  Figure  15.  The  grid  on 
the  steel  surface  was  illuminated  by  an  explosive  argon 
candle  not  shown  in  the  figure.  The  experiments  did  not  pro- 
vide any  velocity  data  due  to  the  low  velocity  of  the  surface. 
This  was  true  even  though  the  grid  pitch  was  twice  that  used 
in  the  feasibility  study.  Later  measurements  using  pins 
(Section  3. 2.3.3)  showed  that  the  velocity  was  only  8 percent 
of  the  value  obtained  with  the  thin  plate. 

Although  the  Moire  fringe  technique  failed  in  these 
particular  experiments,  it  is  a promising  technique  for 
higher  velocities  and  can  be  used  to  measure  velocities  over 
an  entire  surface. 

3. 2. 3. 3 Velocity  Pin  Measurements 

In  the  final  ten  experiments  the  surface  velocity  was 
measured  by  means  of  charged  pins  which  were  shorted  by 
contact  with  the  steel  plate.  These  pins,  of  soft  copper, 
were  held  in  a plexiglas  retainer  and  supported  just  above 
the  plate.  The  velocity  was  calculated  by  dividing  the  known 
spacing  of  the  pins  by  the  time  lapse  between  shorting  signals, 
In  some  cases,  the  velocity  was  not  constant  as  the  plate 
rose,  but  rather  decreased  with  time,  In  these  cases  the 
first  measured  velocity  is  reported. 
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3.3  RESULTS 
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The  data  derived  from  these  experiments  are  tabulate^/ 
in  Table  IV.  Figure  16  shows  data  for  a typical  H configura- 
tion and  Figure  17  shows  data  for  a typical  V configuration. 

A telephoto  picture  of  an  experiment  in  progress  is  shewn 
in  Figure  18. 
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3.4  DISCUSSION 


In  this  section  we  discuss  the  results  of  the  experi- 
ments themselves  and  in  relation  to  the  work  of  Wenzel  and 
lisparzaJ^  Comparison  of  the  experimental  results  with  the 
computational  model  is  deferred. 

3.4.1  effects  of  Nonplanar  Stress  Waves 

The  manganin  foil  stress  gauge  is  designed  to  measure 
stress  in  a planar  wave  traveling  normal  to  the  gauge  plane. 

The  waves  studied  in  this  work  do  not  perfectly  meet  the  re- 
quirements of  planarity  and  normality.  This  is  especially  true 
for  gauges  used  to  measure  stress  off  the  axis  of  the  charge. 
However,  the  effects  of  these  disparities  are  probably  rather 
small  in  the  present  applications.  Figure  19  is  an  example 
of  the  worst  effects  of  nonplanarity.  The  oscillations  prior 
to  arrival  of  the  main  wave  are  due  to  the  diverging  nature 
of  the  wave.  The  records  shown  ir.  Figures  16  and  17  are 
typical  of  almost  all  of  the  data. 

3.4.2  Reproducibility 
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One  of  the  principal  goals  of  the  present  study  was 


2092 

Dry 

1.36 

Sphere 

0 

14.1 

23 

- - 

16 

8.62 

14. 1 

3.1 

>20 

2093 

D 

1.36 

Sphere 

4.31 

14.1 

9.8 

. - 

80 

8.62 

14.1 

4.0 

60 

8.62 

14.1 

3.5 

-- 

2106 

D 

1.36 

Sphere 

0 

14.1 

19 

13 

45 

4.31 

14.1 

8.3 

80 

8.62 

14.1 

4.1 

6S 

8.62 

14.1 

2.2 

55 

211S 

Wet 

2.12 

Sphe re 

0 

14.1 

37 

34 

15(29) 

0 

16.6 

24 

13(26) 

0 

19.2 

16 

17(28) 

0 

21.7 

9.3 

14(25) 

2132 

W 

2.12 

Sphere 

0 

14.1 

35 

18 

12(24) 

4.31 

14.1 

10.6 

13 

2143 

W 

2.12 

Sphere 

0 

14.1 

36 

15 

17(32) 

4.31 

14.1 

21 

13(27) 

8.62 

14.1 

8 

SO 

8.62 

14.1 

4 

46 

2138 

W 

2.12 

Pancake 

0 

14.9 

62 

68 

17 

4.31 

14.9 

36 

13(32) 

2139 

W 

2.12 

Pancake 

0 

14.9 

76 

65 

16 

0 

17.4 

52 

17(20) 

0 

20.0 

36 

20(24) 

0 

22.5 

27 

12(16) 

2144 

W 

2.12 

Pancake 

0 

17.4 

45-6^ 

42 

16(20) 

4.31 

17.4 

29-42 

16(20) 

8.62 

17.4 

11-16 

16 

8.62 

17.4 

11-16 

10(14) 

2140 

D 

1.36 

Pancake 

0 

14.9 

69 

21 

17 

8,62 

14.9 

8.8 

21 

21S0 

D 

1.36 

Pancake 

0 

14.9 

68 

45 

14 

4.31 

14.9 

37 

>9 

8.62 

14.9 

10 

34 

8.62 

14.9 

9.0 

>12 

21S3 

D 

1.36 

Pancake 

0 

14.9 

64 

50 

18 

0 

17.4 

38 

17 

0 

20.0 

25 

22 

0 

22.5 

19 

18 

"Horizontal  distance  from  charge  center. 
^Vertical  distance  from  charge  center. 


c 

d 


Time  duration  to  first  stress  minimum.  Value  in  parentheses 
()  includes  any  "precursor"  or  ramp  (see  Figure  16,  Scope  B2). 

♦ indicates  value  projected  to  2.Scm  from  lower  surface  of 
plate. 
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17(a)  Manganin  gauge  record  at  r = 0,  R = 14.'.)  cm;  Peak  Stress 
(>4  Khar.  Horizontal  Scale  10  ;.sce/div.  Vertical  Scale 
17  kbar/div. 


17(h)  Manganin  gauge  record  at  r = 0,  R = 17.4  cm;  Peak  Stress 
38  kbar.  Horizontal  Scale  do  csec/div.  Vertical  Scale 
'v  7 kbar/div. 


Figure  17.  Representative  experimental  data  from  a V con- 
figuration shot  \sith  a pancake  soil  in  dry  sand, 


17(C)  Mangan  in  jiau^c  record  at  r = 0,  R = 20.0  cm;  Peak  Stress 
2.i  kbar.  Horizontal  Scale  20  ;i  sec/d  i v . Vertical  Scale 
■ 7 kbar/Jiv. 


17(d)  Siangan i n gauge  record  at  r = 0,  R = 22.  S cm;  Peak  Stress 
2S  kbar.  Horizontal  Scale  20  uscc/iliv.  Vertical  Scale 
'*>  0 kbar/div. 
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17(e)  Velocity  pin  record  with  free-surface  velocity  of  a,  50m/sec.  ■ 

Pin  spacing  o.02f>  cm;  Horizontal  Scale  10  asec/di v . ■ 
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I-iyure  19.  lypicul  "worst11  result  showing  effects  of 
nonplana  r have  fronts  on  manganin  gauge 
records.  These  records  were  taken  with 
r/R  = 0.01,  using  a spherical  charge  buried 
under  one  inch  of  wet  soil. 


careful  soil  preparation  techniques.  This  proved  to  be  some- 
what difficult  for  experiments  with  spherical  charges  buried 
in  dry  soil,  primarily  because  of  the  difficulty  of  preparing 
dry  soil  packed  around  spheres.  Even  in  this  case,  peak 
stresses  at  similar  locations  were  within  4 kbar  of  each 
other.  Reproducibility  in  other  cases  was  excellent. 

Figure  20  shows  the  on-axis  stress  histories  for  three 
different  shots  with  a pancake  buried  under  dry  soil.  Not 
only  are  the  peak  stresses  within  ± 4 percent,  but  the  shapes 
of  the  curves  arc  almost  indistinguishable  except  for  the 
high  frequency  noise  probably  caused  by  electrical  pickup. 

3.4.3  Effect  of  Vertical  Pos:;  tion  in  Steel  Plate 

We  expect  the  peak  stress  to  decrease  as  pulses  propa- 
gate up  through  the  plate  for  two  reasons.  First,  as  the  wave 
propagates  it  will  diverge  over  a larger  area.  This  effect 
should  be  expected  to  produce  a decrease  in  peak  stress  inversely 
proportional  to  the  souare  of  the  distance  from  the  charge  center 
for  a spherical  geometry.  Second,  since  the  relief  wave  following 
the  peak  stress  will  travel  faster  than  the  loading  wave, 
the  peak  should  decrease  in  magnitude  even  in  the  absence  of 
divergence.  This  effect  of  finite  amplitude  waves  should  be 
more  pronounced  at  higher  stresses. 

The  three  V- configuration  experiments  were  undertaken 
to  quantify  the  expectations  just  expressed.  The  results  are 
shown  graphically  in  Figures  21  and  22.  In  Figure  21,  the 
decay  with  propagation  is  evident.  From  Figure  22  we  can  see 
that  the  decay  goes  approximately  as  the  inverse  square  of 
the  radius.  It  appears  that  the  second  mechanism  of  attenua- 
tion is  unimportant  for  both  spherical  charges  and  pancake  charges. 
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Although  we  did  not  measure  impulse  directly  in  our 
experiments  it  can  be  determined  from  the  area  under  a pressure- 
time curve.  We  note  that  the  duration  and  shape  of  the  main 
pressure  pulse  does  not  change  significantly  as  the  pulse 
propagates;  however,  the  amplitude  does  decrease.  Consequently, 
the  specific  impulse  decreases  on  propagation  in  the  same 
manner  as  the  peak  stress. 


3.4.4  Dependence  on  Distance  off  Axis 
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The  peak  stress  is  a very  sensitive  function  of  the 
location  of  the  gauge  relative  to  the  symmetry  axis  of  the 
experiment.  This  is  especially  true  for  pancake  charges. 

This  can  be  clearly  seen  in  Figure  23  where  peak  stress  is 
shown  as  a function  of  the  dimensionless  variable  r/R  where 
r and  R are  the  horizontal  and  vertical  distances  from  the 
charge  center  respectively.  This  effect  may  account  for  most 
of  the  scatter  in  the  data.  An  error  in  locating  the  plate 
over  the  charge  axis  of  1 cm  can  cause  a change  in  the  peak 
stress  of  as  much  as  25  percent  of  the  peak  stress  directly 
above  the  charge. 


£ I It  should  be  noted  that  all  of  the  measurements  reported 

? here  were  taken  at  distances  off  axis  of  less  than  the  charge 

\ diameter. 
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3.4.5  Effect  of  Soil  State 


As  we  mentioned  in  Subsection  3.1.2,  Wenzel  and  Esparza ^ 
concluded  that  bu. ial  did  little  to  increase  peak  stress 
from  pancake  charges  (relative  to  charges  detonated  in  air) , 
but  that  the  soil  had  a considerable  effect  for  spherical 
charges.  We  argued  that  the  mechanical  impedance  of  the  soil 
should  be  the  important  criterion  leading  to  any  difference 
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in  peak  stress  for  charges  buried  or  in  air.  If  this  were 
true,  we  should  expect  to  see  a much  greater  peak  stress  for 
spherical  charges  under  our  "wet"  soil  than  under  our  "dry" 
soil.  Indeed,  this  is  the  case,  as  can  be  seen  from  Figure  23 
and  Table  IV.  For  the  pancake  charges,  on  the  other  hand, 
we  find  no  significant  differences  in  peak  stress,  just  as 
we  would  expect  from  the  observations  of  Wenzel  and  Esparza. ^ 
These  results,  therefore,  seem  to  confirm  their  observations 
and  our  explanation  of  them.  This  is  at  variance,  however, 
with  their  observation  that  pancakes  buried  in  mud  produce 
much  higher  peak  stresses  than  those  in  dry  soil.  The  question 
of  why  we  see  no  effect  of  burial  of  pancakes  is  not  easy  to 
answer,  but  if  there  is  little  difference  between  charges 
in  air  and  those  in  dry  soil,  we  should  expect  little  variation 
between  those  in  two  different  soils. 

3.4.6  Comparison  with  Earlier  Work 

The  broad  features  of  the  present  work  agree  quite 
well  with  those  seen  by  Wenzel  and  Esparza^  As  just  men- 
tioned, our  findings  on  the  effects  of  different  soils  agree 
with  theirs  on  soils  versus  air.  Our  study  confirms  the 
importance  of  soil  impedance  for  spheres,  as  inferred  by 
their  modeling  relations. 

However,  there  is  one  glaring  difference  between  their 
results  and  ours.  That  is  the  magnitude  of  the  stresses 
involved.  Stresses  observed  in  the  present  study  were  uni- 
formly about  three  times  greater  than  those  in  corresponding 
experiments  by  Wenzel  and  Esparza.  ^ The  explosives  used 
(C-4  in  the  present  work,  pentolite  in  the  other^  are  not 
sufficiently  different  to  account  for  this  effect.  The 
pulses  seen  in  the  earlier  study  were  of  much  longer  duration 
than  those  we  saw.  If  we  estimate  the  specific  impulse  at 
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the  lower  surface  of  the  plate  from  our  pressure- time  data, 
we  find  that  impulses  are  fairly  close  in  the  two  sets  of 
experiments. 

We  tentatively  ascribe  the  differences  to  a combination 
of  geometric  dispersion  and  finite  wave  amplitude  effects  in 
the  bar  gauges.  Note  from  our  V configuration  results  of 
Figure  21  that  our  measurements  at  4- inch  depths  in  the  plates 
give  peak  stresses  close  to  the  SwRI  results.  The  different 
geometries  of  the  plates  and  bars  prevent  any  quantitative 
comparison  of  waveforms  at  a given  distance  from  the  impact 
surfaces.  The  bar  gauges  do  not  give  a stress  history  that 
can  be  related  to  a particular  depth  in  a plate.  The  1- inch- 
deep  manganin  gauge  results  should  be  closer  to  impact  stress. 

3.5  SUMMARY  AND  CONCLUSIONS 


Twelve  experiments  were  carried  out  with  spherical 
and  pancake  charges  buried  under  wet  and  dry  soils.  In  nine 
of  these  experiments  stress  histories  were  measured  in  a plane 
within  a steel  plate  suspended  over  the  charge.  One  of  these 
gauges  was  directly  over  the  charge  and  the  others  were  located 
at  various  radii  from  the  perpendicular  axis  through  the  plate 
from  the  charge  center.  In  three  of  the  experiments  the  stress 
histories  were  measured  at  four  locations  along  the  axis  within 
the  plate.  Careful  attention  to  soil  preparation  led  to  very 
good  reproducibility  of  results. 
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The  experimental  investigations  led  to  the  following 
conclusions : 


• The  blast  effect  from  a buried  sphere  of 
high  explosive  depends  on  the  soil 
condition. 


The  blast  effect  from  a buried  pancake  of 
high  explosive  is  independent  of  the  soil 
condition . 


• The  peak  stress  and  specific  impulse  fall 
off  approximately  as  1/R2  in  a steel  plate 
over  a buried  charge. 


The  magnitude  of  peak  stresses  measured  in 
the  present  study  is  about  three  times 
greater  than  those  measured  previously  at 
SwRI . However,  the  impulses  in  the  two 
studies  are  similar. 
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IV.  CALCULATIONS 

Three  buried  charge  calculations  were  performed  using 
a slightly  modified  version  of  the  numerical  model  previously 
reported. 1 * 1 The  additional  modifications  are  reported  in 
Appendix  A.  Calculation  1 involved  a spherical  charge  and 
was  modelled  after  Experiments  2115,  2132  and  2143  of  Table  IV. 
Calculation  2 involved  a disc-shaped  or  pancake  charge  and 
was  modelled  after  Experiments  2138,  2139  and  2144  of  Table  IV. 
Calculation  3 was  not  modelled  after  any  of  the  experiments 
and  involved  a steel  plate  in  contact  with  the  soil  surface. 

The  results  of  these  three  calculations  are  reported  here 
and  comparisons  with  experiment  are  made  whenever  possible. 

For  more  details  regarding  the  numerical  model,  Refs.  1 and  2 
and  Appendix  A of  this  report  should  be  consulted. 

Figure  24  shows  a typical  initial  configuration  employed 
in  the  buried  charge  calculations.  The  values  of  depth-of- 
burial,  n,  and  standoff,  S,  employed  in  the  three  calculations 
are  given  in  Table  V.  Tracer  particles,  referred  to  as 
sensors  in  Figure  24,  were  placed  at  various  positions  through- 
out the  steel  plate  for  the  purpose  of  monitoring  the  average 
axial  normal  stress  components  at  those  locations  as  a function 
of  time  and  comparing  the  predicted  results  with  experiment. 

In  the  calculations,  rows  of  tracer  particle  sensors  were 
placed  along  the  front  surface  of  the  steel  plate  and  at 
one- inch  depth  intervals,  in  the  experimental  set-up,  gauges 
for  measuring  the  stress  were  embedded  in  the  plate  at  similar 
positions;  however,  no  gauges  were  placed  on  the  front  of  the 
plate  since  they  would  have  been  damaged  by  the  impacting  debris. 


Figure  25  shows  the  predicted  configuration  of  the  high 
explosive  charge,  the  soil  and  the  steel  plate  at  various 
times  after  initiation  for  a typical  solution  of  the  type 
that  is  of  interest  here.  Since  the  configurational  plots 
are  quite  similar  and  provide  only  qualitative  information, 
they  will  not  be  repeated  for  each  of  the  three  solutions. 
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I’ig.  25- -Configurations  at  various  times  after  initiation 
for  a typical  buried  mine-steel  plate  interaction 
calculation. 
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The  predicted  total  normal  axial  stress  pulse  as  a 
function  of  distance  into  the  plate  is  shown  at  various  times 
in  Figure  26  for  Calculation  1 of  Table  V.  Initiation  of 
the  buried  charge  occurred  at  t = 0. 

Figure  27  is  a plot  of  the  predicted  axial  stress 
history,  o^Ct),  at  the  front  surface  of  the  steel  plate 
and  at  various  gage  locations  within  the  plate.  It  is 
seen  in  Figure  27  that  the  stress  pulse  attenuated  quite 
rapidly  as  it  propagated  into  the  plate.  The  symbols,  (x) , 
in  Figure  27  show  the  experimentally  measured  peak  stresses 
at  the  various  gage  locations.  It  is  seen  that  agreement 
between  theory  and  experiment  is  quite  good. 

In  Figure  28  the  predicted  axial  stress  histories, 
Oyy(t),  at  gages  located  one  inch  into  the  plate  and  at 
various  radial  distances,  r,  are  shown  for  Calculation  1. 

The  experimental  data  points  given  by  the  symbol  (x)  indicate 
the  measured  peak  stresses  at  the  various  gage  locations. 

As  is  seen  from  Figure  28,  while  there  is  some  spread  in 
the  experimental  data,  agreement  between  theory  and  experi- 
ment is  quite  good. 

Stress  plots  similar  to  those  shown  in  Figures  26 
through  28  are  provided  in  Figures  29  through  31  for 
Calculation  2 of  Table  vr.  In  addition  to  those  plots,  stress 
histories  at  gages  located  two  inches  deep  in  the  steel 
plate  are  provided  for  Calculation  2 in  Figure  32.  As  seen 
in  Figure  30  the  predicted  stress  attenuation  through  the 
thickness  of  the  plate  is  in  good  agreement  with  the 
measured  stress  attenuation;  however,  the  predicted  stress 
peaks  are  consistently  lower  than  the  measured  values  by 
a few  per  cent. 
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I-'igure  28. 


Plot  of  predicted  axial  stress  history,  Oyy(t), 
at  gauges  located  one  inch  into  the  plate  and  at 
various  radial  distances,  r,  for  Calculation  1. 
The  experimentally  measured  values  are  shown  by 
the  symbol  (x) . 
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The  predicted  axial  stress  attenuation  in  the  radial 
direction  for  Calculation  2 does  not  agree  as  well  with 
experiment  as  does  the  attenuation  through  the  plate.  As 
seen  in  Figures  31  and  32,  the  predicted  radial  attenuation 
is  not  as  great  as  the  measured  values.  The  reason  for 
this  difference  in  the  theoretical  and  experimental  curves 
for  the  radial  attenuation  of  the  peak  axial  stress  components 
for  Calculation  2 is  not  completely  understood  at  this  time; 
however,  it  is  apparent  from  the  last  plot  in  Figure  3 that 
the  predicted  debris  momentum  did  not  decrease  monotonically 
in  the  radial  direction. 

Calculation  3 involved  a steel  plate  in  contact  with 
the  ground  and  the  predicted  total  normal  axial  stress  as 
a function  of  distance  into  the  plate  is  plotted  at  various 
times  after  charge  initiation  in  Figure  33.  Figure  34 
shows  the  predicted  attenuation  of  the  peak  stress  as  the 
wave  propagates  through  the  plate.  The  experimental  matrix 
did  not  include  a configuration  similar  to  Calculation  3 
so  that  no  data  are  currently  available  for  comparison 
with  the  predicted  peak  stresses.  Figures  35  and  36  show 
the  attenuation  of  the  axial  stress  wave  as  it  propagates 
into  the  plate  at  radial  distances  of  1.8  and  3.7  cm,  respec- 
tively, for  Calculation  3.  While  no  experimental  data  were 
available  for  comparison  with  the  results  of  Calculation  3, 
these  results  were,  however,  compared  with  results  obtained 
at  MI:.RDC 1 2 1 3 jn  which  the  Lagrangian  code  TGODY^^  was 


employed.  The  predicted  values  of  the  impulse  delivered  to 

the  steel  plate  along  the  axis  of  symmetry  differed  by  only 

3 [21] 

The  predicted  total  momentum  given  to  the  steel  plate 
as  a function  of  time  after  debris  impact  is  plotted  in 
Figure  37  for  each  of  the  three  calculations  of  Table  V. 


Interface 


5 Numerical  techniques  have  been  developed  and  proven 

successful  for  predicting  the  effects  of  buried  mines  on 
simple  structures  located  above  the  ground,  lixperimental 
techniques  for  determining  dynamic  material  properties  of 
soils  have  been  employed  to  obtain  equation  of  state  para- 
meters which  are  necessary  as  part  of  the  input  to  the  numeri- 
cal model . 


Three  two-dimensional  calculations  were  performed. 

In  these  calculations  the  buried  charges  were  detonated, 
the  expanding  detonated  products  accelerated  the  soil  and 
the  debris  impacted  a steel  plate  located  above  ground 
level.  The  stress  distributions  throughout  the  steel  plates 
as  a function  of  time  after  charge  initiation  were  determined 
and  compared  with  data  from  actual  tests.  The  agreement 
between  theory  and  experiment  was  quite  good. 

Specific  conclusions  regarding  the  results  of  the 
dynamic  soil  properties  tests  and  the  actual  buried  charge- 
steel  plate  interaction  tests  are  listed  at  the  ends  of 
Sections  11  and  III,  respectively. 

In  general  it  can  be  said  that  the  combined  theoretical - 
experimental  program  has  significantly  advanced  the  state-of- 
the-art  of  predictive  techniques  related  to  the  effects  of 
buried  charges  on  simple,  above-ground  structures. 
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APPENDIX  A 

HELP  CODE  MODIFICATIONS 

Some  changes  have  been  made  to  the  version  of  the  HELP 
code  reported  in  Reference  2.  The  major  change  has  been  to 
incorporate  an  improved  artificial  viscosity  formulation. 

The  artificial  viscosity  which  is  added  to  the  stress  at  the 
boundaries  between  neighboring  cells  K and  K1  is  of  the  form 
Q=QLIN*^*U"AV,  where  QLIN  is  an  input  parameter,  p = 0.5  (p(K) 

+p  (Kl))  where  p is  the  cell  density,  C = 0.5  (C(K)  + C (K 1 ) ) 
where  C is  the  cell  sound  speed,  and  AV  is  the  velocity 
difference  between  the  two  cells.  For  the  top  boundary  the 
velocity  difference  is  AV  = VfX;  - V(K1)  where  V is  the  cell 
axial  velocity,  while  for  the  right  boundary  it  is  AV  = U(K)  - 
U (Kl)  where  U is  the  cell  radial  velocity. 

This  formulation  has  required  the  addition  of  a func- 
tion, SNDSPD , which  calculates  the  cell  sound  speed.  For  the 
high  explosive  package  the  sound  speed  is  gi/en  by  c=/v  (y- 1 ) t'. , 
where  E is  the  specific  internal  energy  and  y is  an  input 
parameter.  For  other  materials  the  sound  speed  is  giver,  by 

C = C -*0.5  /P,  where  P is  the  cell  pressure  and  C is 
0 0 
the  bulk  sound  speed.  For  mixed  cells,  a volume  weighted 

average  is  used. 

Two  input  ariables  have  been  adue.  Lo  the  Z block. 

The  first  of  these,  in  location  11,  is  QL1,  the  artificial 
viscosity  coefficient.  The  other  is  an  integer,  NVRTE,  in 
location  107.  This  variable  must  be  defined  eqyal  to  NVRTEX 
in  the  initial  problem  setup. 

Another  code  modification  was  made  in  order  to  reduce 
the  zoning  requirements  in  che  vicinity  of  the  debris-steel 
plate  interface.  in  the  unmodified  version  of  the  code, 
the  under-dense  soil  entering  the  mixed,  interface  cell  dur- 
ing the  transport  phase  transferred  momentum  to  the  plate 
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for  several  cycles  before  a significant  stress  gradient  acted 
across  the  interface.  The  result  of  this  lack  of  synchroni- 
zation was  to  artifically  broaden  the  stress  pulse  delivered 
to  the  plate.  This  problem  could  have  been  eliminated  by 
choosing  finer  zoning;  however,  the  cost  of  the  calculations 
would  then  have  been  significantly  greater.  The  modification 
that  was  made  simply  inhibited  the  transport  into  the  steel 
cells  until  a significant  stress  gradient  developed  across 
the  interface  so  that  the  Phase  II  (mass  transport)  and  Phase 
III  (stress  gradient)  effects  would  be  synchronized. 
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